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The best particle candidates for non-baryonic cold dark matter are reviewed, namely, neutralino, axion, axino 
and Majoron. These particles are considered in the context of cosmological models with the restrictions given by 
the observed mass spectrum of large scale structures, data on clusters of galaxies, age of the Universe etc. 



1. Introduction (Cosmological environ- 
ment) 

Presence of dark matter (DM) in the Universe 
is reliably established. Rotation curves in many 
galaxies provide evidence for large halos filled by 
nonluminous matter. The galaxy velocity distri- 
bution in clusters also show the presence of DM in 
intercluster space. IRAS and POTENT demon- 
strate the presence of DM on the largest scale in 
the Universe. 

The matter density in the Universe p is usually 
parametrized in terms of f2 = p/pc, where pc ~ 
1.88 • 10"^^ g/cm? is the critical density and 
h is the dimensionless Hubble constant defined 
es h ~ Hq/ {100km. s^^.Mpc~^). Different mea- 
surements suggest generally 0.4 < h < 1. The 
recent measurements of extragalactic Cepheids in 
Virgo and Coma clusters narrowed this interval to 
0.6 < ft. < 0.9. However, one should be cautious 
about the accuracy of this interval due to uncer- 
tainties involved in these difficult measuremets. 

Dark Matter can be subdivided in baryonic 
DM, hot DM (HDM) and cold DM (CDM). 

The density of baryonic matter found from nu- 
cleosynthesis is given § as Qth^ = 0.025 ±0.005. 

Hot and cold DM are distinguished by velocity 
of particles at the moment when horizon crosses 
the galactic scale. If particles are relativistic they 
are called HDM particles, if not - CDM. The nat- 
ural candidate for HDM is the heaviest neutrino, 



most naturally r-neutrino. Many new particles 
were suggested as CDM candidates. 

The structure formation in Universe put strong 
restrictions to the properties of DM in Universe. 
Universe with HDM plus baryonic DM has a 
wrong prediction for the spectrum of fluctuations 
as compared with measurements of COBE, IRAS 
and CfA. CDM plus baryonic matter can ex- 
plain the spectrum of fluctuations if total density 
Qa ~ 0.3. 

There is one more form of energy density in the 
Universe, namely the vacuum energy described 
by the cosmological constant A. The correspond- 
ing energy density is given by flA — A/{3Hq). 
Quasar lensing and the COBE results restrict the 
vacuum energy density: in terms of i}\ it is less 
than 0.7. 

Contribution of galactic halos to the total den- 
sity is estimated as il ~ 0.03 — 0.1 and clusters 
give ft « 0.3. Inspired mostly by theoretical mo- 
tivation (horizon problem, flatness problem and 
the beauty of the inflationary scenarios) fio = 1 
is usually assumed. This value is supported by 
IRAS data and POTENT analysis. No observa- 
tional data significantly contradict this value. 

There are several cosmological models based on 
the four types of DM described above (baryonic 
DM, HDM, CDM and vacuum energy). These 
models predict different spectra of fluctuations to 
be compared with data of COBE, IRAS, CfA etc. 
They also produce different effects for cluster- 



2 



cluster correlations, velocity dispersion etc. The 
simplest and most attractive model for a cor- 
rect description of all these phenomena is the 
so-called mixed model or cold-hot dark matter 
model (CHDM). This model is characterized by 
following parameters: 

= 0,^0 = + ^CDM + ^HDM = 1, 

Hq « 50 kms-^Mpc-^{h « 0.5), 
^CDM ■■ ^HDM ■■ ~ 0.75 : 0.20 : 0.05, (1) 

where VIhdm ~ 0.2 is obtained in ref. 1^ from 
damped Lya data. Thus in the CHDM model 
the central value for the CDM density is given by 

^CDMh" - 0.19 (2) 

with uncertaities within 0.1. 

The best candidate for the HDM particle is r- 
neutrino. In the CHDM model with = 0.2 
mass of r neutrino is m^^ k, 4.7 eV . This com- 
ponent will not be discussed further. 

The most plausible candidate for the CDM par- 
ticle is probably the neutralino (x): it is massive, 
stable (when the neutralino is the lightest super- 
symmeric particle and if R-parity is conserved) 
and the xX'^'Unihilation cross-section results in 
D.^h?' - 0.2 in large areas of the neutralino pa- 
rameter space. 

In the light of recent measurements of the Hub- 
ble constant the CHDM model faces the age prob- 
lem. The lower limit on the age of Universe 
io > 13 Gyr (age of globular clusters) imposes the 
upper limit on the Hubble constant in the CHDM 
model < 50 kms~^Mpc~^. This value is in 
slight contradiction with the recent observations 
of extragalactic Cepheids, which can be summa- 
rized as Hq > 60 kms^^Mpc~^. However, it is 
too early to speak about a serious conflict tak- 
ing into account the many uncertainties and the 
physical possibilities (e.g. the Universe can be 
locally overdense - see the discussion in ref . Q ) . 

The age problem, if to take it seriously, can 
be solved with help of another successful cosmo- 
logical model ACDM. This model assumes that 
fio = 1 is provided by the vacuum energy (cos- 
mological constant A) and CDM. From the limit 
Ha < 0.7 and the age of Universe one obtains 
^CDAi > 0.3 and h < 0.7. Thus this model 



also predicts ^IcDAih^ « 0.15 with uncertainties 
0.1. Finally, we shall mention that the CDM with 
rio = ^CDM = 0.3 and h = 0.8, which fits the ob- 
servational data, also gives flh^ « 0.2. Therefore 
O/i^ « 0.2 can be considered as the value common 
for most models. 

In this paper we shall analyze several candi- 
dates for CDM, best motivated from the point of 
view of elementary particle physics. The motiva- 
tions are briefly described below. 
Neutralino is a natural lightest supersymmetric 
particle (LSP) in SUSY. It is stable if R-parity is 
conserved, il^ ^ ^cdm is naturally provided by 
annihilation cross-section in large areas of neu- 
tralino parameter space. 

Axion gives the best known solution for strong 
CP-violation. Cla ^ ^cdm for natural values of 
parameters. 

Axino is a supersymmetric partner of axion. It 
can be LSP. 

Majoron is a Goldstone particle in spontaneously 
broken global U{1)b-l or U{1)l- KeV mass can 
be naturally produced by gravitational interac- 
tion. 

Apart from cosmological acceptance of DM 
particles, there can be observational confirma- 
tion of their existence. The DM particles can 
be searched for in the direct and indirect experi- 
ments. The direct search implies the interaction 
of DM particles occurring inside appropriate de- 
tectors. Indirect search is based on detection of 
the secondary particles produced by DM particles 
in our Galaxy or outside. As examples we can 
mention production of antiprotons and positrons 
in our Galaxy and high energy gamma and neu- 
trino radiation due to annihilation of DM parti- 
cles or due to their decays. 

2. Axion 

The axion is generically a light pseudoscalar 
particle which gives natural and beautiful solu- 
tion to the CP violation in the strong interaction 
Q (for a review and references seeQ). Sponta- 
neous breaking of the PQ-symmetry due to VEV 
of the scalar field < 4> >= fpg results in the 
production of massless Goldstone boson. Though 
fpQ is a free parameter, in practical applications 
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it is assumed to be large, /pq ~ 10^" — 10^^ QeV 
and therefore the PQ-phase transition occurs in 
very early Universe. At low temperature T ^ 
^QCD ~ 0.1 GeV the chiral anomaly of QCD 
induces the mass of the Goldstone boson ma ^ 
^qcd/ fPQ ■ This massive Goldstone particle is 
the axion. The interaction of axion is basically 
determined by the Yukawa interactions of field(s) 
(j) with fermions. Triangular anomaly, which pro- 
vides the axion mass, results in the coupling of 
the axion with two photons. Thus, the basic for 
cosmology and astrophysics axion interactions are 
those with nucleons, electrons and photons. 
Numerically, axion mass is given by 

ma = 1.9 • 10-3(7V/3)(10i° GeV/fpg) eV, (3) 

where A'^ is a color anomaly (number of quark 
doublets). 

All coupling constants of the axion are inversely 
proportional to fpg and thus are determined by 
the axion mass. Therefore, the upper limits on 
emission of axions by stars result in upper limits 
for the axion mass. Of course, axion fluxes can- 
not be detected directly, but they produce addi- 
tional cooling which is limited by some observa- 
tions (e.g., age of a star, duration of neutrino pulse 
in case of SN etc). In Table 1 we cite the upper 
limits on axion mass from ref.||], compared with 
revised limits, given recently by Raffelt |6) . 



Table 1 

Astrophysical upper limits on axion mass 





1990 |] 


1995 § 


sun 


1 eV 


1 eV 


red giants 


1 ■ 10-2 eV 


very 

uncertain 


hor .-branch 


not 


0.4 eV 


stars 


considered 


SN 1987A 


1 • 10-^ eV 


1 • 10-2 eV 



As one can see from the Table the strong upper 
limit, given in 1990 from red giants, is replaced 
by the weaker limit due to the horizontal-branch 



stars. The upper limit from SN 1987A was re- 
considered taking into account the nucleon spin 
fluctuation inA^-fA^ — > A^4-A^-|-a axion emission. 

There are three known mechanisms of cosmo- 
logical production of axions. They are (i)thermal 
production, (i) misalignment production and (iii) 
radiation from axionic strings. 

The relic density of thermally produced axions 
is about the same as for light neutrinos and thus 
for the mass of axion ma ~ 10-2 eV this compo- 
nent is not important as DM. 

The misalignment production is clearly ex- 
plained in ref.|^. 

At very low temperature T <^ ^qcd the mas- 
sive axion provides the minimum of the potential 
at value 9 = 0, which corresponds to conservation 
of CP. At very high temperatures T :s> A^qcd 
the axion is massless and the potential does not 
depend on 9. At these temperatures there is no 
reason for 9 to be zero: its values are different in 
various casually disconnected regions of the Uni- 
verse. When T — > Aqcd the system tends to go 
to potential minimum (at 9^0) and as a re- 
sult oscillates around this position. The energy 
of these coherent oscillations is the axion energy 
density in the Universe. From cosmological point 
of view axions in this regime are equivalent to 
CDM. The energy density of this component is 
approximately 

Qah^ « 2 • (ma/10-5 el/)-i i*. (4) 

Uncertainties of the calculations can be estimated 
as 10±0-5. 

Axions can be also produced by radiation of 
axionic strings |^,|8|. Axionic string is a one- 
dimension vacuum defect < <ppQ >= 0, i.e. a 
line of old vacuum embedded into the new one. 
The string network includes the long strings and 
closed loops which radiate axions due to oscilla- 
tion. There were many uncertainties in the axion 
radiation by axionic strings (see ref. |5| for a re- 
view). Recently more detailed and accurate cal- 
culations were performed by Battye and Schellard 
1^ . They obtained for the density of axions 

riah^ « A(m JlQ-^ eV)-^-^^ (5) 

with A limited between 2.7 and 15.2 and with 
uncertainties of the order 10^°'^. The overpro- 
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duction condition flah^ > 1 imposes lower limit 
on axion mass ma > 2.3 • 10~^ eV. 

Fig.l shows the density of axions ^lah^ as a 
function of the axion mass ma- The upper limits 
on axion mass from Table I are shown above the 
upper absciss (limits of 1990) and below lower ab- 
sciss (limits of 1995). The overproduction region 
V,ah^ > 1 and the regions excluded by astrophysi- 
cal observations |^ are shown as the dotted areas. 

The axion window of 1995 (shown as undot- 
ted region) became wider and moved to the right 
as compared with window 1990. The horizontal 
strip shows ^Icdm — 0.2 ±0.1 as it was discussed 
in Introduction. One can see from Fig.l that 
string and misaligment mechanisms provide the 
axion density as required by cosmological CDM 
model, if axion mass is limited between 7-10~^ eV 
and 7 • 10"'* eV. However, in the light of un- 
certainties, mostly in the calculations of axion 
production, one can expect that this " best calcu- 
lated" window is between 3 • 10~^ and 10"'^ eV. 
This region is partly overlapped with a possible 
direct search for the axion in nearest-future ex- 
periments (see Fig.l and refs.Q). 

3. Axino 

In supersymmetric theory the PQ-solution 
for strong CP-violation should be generalized. 
Within this theory the PQ symmetry breaking re- 
sults in the production of the Goldstone chiral su- 
permultiplet which contains two scalar fields and 
their fermionic partner - axino (a). The scalar 
fields enter the supermultiplet in the combination 
ifpQ + s) exp(a//pQ), where s is a scalar field, 
saxino, which describes the oscillations of the ini- 
tial field (j) around its VEV value < (j) >= fpg, 
and a is the ax;ion field. This phase transition in 
the Universe occurs at temperature T ~ fpQ- As 
we saw in the previous section the axion is mass- 
less at this temperature and since supersymme- 
try is not broken yet, the axino and saxino are 
massless, too. The axion acquires the mass in the 
usual way due to chiral anomaly at T Aqcd, 
while saxino and axino obtain the masses due to 
global supersymmetry breaking. 

The saxino is not of great interest for cosmol- 
ogy: it is heavy and it decays fast (mostly into 



two gluons). The axino can be the lightest super- 
symmetric particle and thus another candidate 
for DM. 

How heavy the axino can be? The mass of ax- 
ino has a very model dependent value. In the phe- 
nomenological approach, using the global super- 
symmetry breaking parameter Msusy one typi- 
cally obtains (e.g. 0,|l|) 

™a ^ Msusy/ fpQ (6) 

For example, if global SUSY breaking occurs due 
to VEV of auxiliary field of the goldstino su- 
permultiplet < F >= Fg, then the axino mass 
appears due to interaction term [g/ fpQ)aaF (F 
has a dimension M^), and using < F >= Fg = 
Msusy arrives at the value (||). 

The situation is different in supergravity. In 
rcf . Il^ the general analysis of the axino mass is 
given in the framework of local supersymmetry. 
It was found that generically the mass of axino in 
these theories is rria ^ ~ 100 GeV. Even in 
case when axino mass is small at tree level, the 
radiative corrections raise this mass to the value 
^ 'm-3/2- This result holds for the most general 
form of superpotential.The global SUSY result, 
'rria ^ '^3/2//pQi t)® reproduced in the local 
SUSY only if one of the superpotential couphng 
constants is very small, A < 10~^, which implies 
fine-tuning. Thus, the axino is too heavy to be a 
CDM particle. 

The only exceptional case was found by Goto 
and Yamaguchi |Q. They demonstrated that in 
case of no-scale superpotential the axino mass 
vanishes and the radiative corrections in some 
specific models can result in the axino mass 
10 — 100 keV, cosmologically interesting. This 
beautiful case gives essentially the main founda- 
tion for axino as CDM particle. 

The cosmological production of axinos can oc- 
cur through thermal production fl^ or due to 
decays of the neutralinos , ||l^ . The axion chi- 
ral supermultiplet contains two particles which 
can be CDM particles, namely axion and ax- 
ino. In this section we are interested in the 
case when axino gives the dominant contribu- 
tion. In particular this can take place in the range 
2 • 10^ GeV < fpQ < 2.7 • 10^" GeV where axions 
are cosmologically unimportant. 
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Figure 1. Axion window 1995. The curves "therm." and "misalign." describe the thermal and misaligne- 
ment production of axions, respectively. The dash-dotted curve corresponds to the calculations by Davis 
10 1 for string production. The recent refined calculations are shown by two dashed lines for two 
extreme cases, respectively. The other explanations are given in the text. 



Since axino interacts with matter very weakly, 
the decoupling temperature for the thermal pro- 
duction is very high |l6| : 

Trf « 10^ Gey(/PQ/10" GeV). (7) 

Therefore, axinos are produced thermally at the 
reheating phase after inflation. The relic concen- 
tration of axinos can be easily evaluated for the 
reheating temperature Tr as 

o ^2 ^ n , rn-a ^ 3 -10^° Gey ,2 Tr 

Reheating temperature Tr < 10^ GeV gives no 
problem with the gravitino production. The relic 
density (||) provides VtcDAih^ ^ 0.2 for a reason- 
able set of parameters ma, fpQ and Tp. One can 
easily incorporate in these calculations the addi- 
tional entropy production if it occurs at EW scale 



If the axino is LSP and the neutralino is the 
second lightest supersymmetric particle, the ax- 
inos can also be produced by neutralino decays 
l|ll],il,0- According to estimates of ref.0 
the axinos are produced due to x — > 5, + 7 de- 
cays at the epoch with red-shift Zdec ^ 10^. Ax- 
inos are produced in these decays as ultrarela- 
tivistic particles and the free-streeming prevents 
the growth of fluctuations on the horizon scale 
and less. At red-shift z„r ~ lO'* axinos become 
non-relativistic due to adiabatic expansion (red 
shift). From this moment on the axinos behave as 
the usual CDM and the fluctuations on the scales 
A > (1 -|- Znr)ctnr (which Correspond to a mass 
larger than IO^^Mq) grow as in the case of stan- 
dard CDM. For smaller scales the fluctuations, 
as was explained above, grow less than in CDM 
model. Therefore, as was observed in ref.||l^, the 
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axinos produced by neutralino decay behave like 
HDM. It means that axinos can provide generi- 
cally both components, CDM and HDM, needed 
for description of observed spectrum of fluctua- 
tions. 

Unfortunatelly stable axino is unobservable. In 
case of very weak R-parity violation, decay of ax- 
inos can produce a diffuse X-ray radiation, with 
practically no signature of the axino. 

4. Majoron 

The Majoron is a Goldstone particle associ- 
ated with spontaneously broken global U{1)b-l 
or U{1)l symmetry. The symmetry breaking oc- 
curs due to VEV of scalar field, < a Vg, and 
a splits into two fields, p and J: 

a ^ {vs + p)e-Kp{iJ/vs). (9) 

The field J is the Majoron. A mass of ^ keV 
can be obtained due to gravitational interac- 
tion ||l|l,[|9). The keV Majoron has the great 
cosmological interest since the Jeans mass associ- 
ated with this particle, m jeans ~ "^pz/^jj gives 
the galactic scale M ~ lO^^M©. In all other re- 
spects the keV Majoron plays the role of a CDM 
particle. It is assumed usually that the Majoron 
interacts directly only with some very heavy par- 
ticles (e.g. with the right-handed neutrino vr). 
It results in very weak interaction of the Majoron 
with the ordinary particles (leptons, quarks etc) 
and thus makes the Majoron "invisible" in the 
accelerator experiments. 

The cosmological production of the Majoron 
occurs through thermal production and ra- 

diation by the strings as in case of the ax- 
ion. However, under imposed observational con- 
straints, the Majoron in models[^,|jl^ has to be 
in general unstable with lifetime much shorter 
than the age of Universe. A successful model 
was developed in ref.||2l[], where the Majoron was 
assumed to interact with the ordinary particles 
through new heavy particles. For the cosmologi- 
cal production it was considered the phase transi- 
tion associated with the global U{1)b-l symme- 
try breaking, when the Majorons were produced 
both directly and through p ^ J + J decays. 



Another interesting possibility was considered 
recently in ref.|^^. In this model the Majoron 
is rather strongly coupled with and Vt neu- 
trinos {Jv^Vt coupling). The Majorons are pro- 
duced through Vr ^ J + Vi_i decays. The strong 
interactions between Majorons reduces the relic 
abundance of the Majorons to the cosmologically 
required value. 

The Majoron signature in all these models is 
given by J 7 -I- 7 decays, which result in the 
production of keV X-ray line in the X-ray back- 
ground radiation. 

5. Neutralino 

The neutralino is a superposition of four spin 
1/2 neutral fields: the wino IV3, bino B and two 
Higgsinos Hi and H2: 

X ^ C1W3 + C2B + C3H1 + C4H2 (10) 

The neutralino is a Majorana particle. With a 
unitary relation between the coefficients Ci the 
parameter space of neutralino states is described 
by three independent parameters, e.g. mass of 
wino M2, mixing parameter of two Higgsinos /i, 
and the ratio of two vacuum expectation values 
tan/3 = V2/V1. 

In literature one can find two extreme ap- 
proaches describing the neutralino as a DM par- 
ticle. 

(i) Phenomenological approach. The allowed 
neutralino parameter space is restricted by the 
LEP and CDF data. In particular these data 
put a lower limit to the neutralino mass, m-^ > 
20 GeV. In this approach only the usual GUT re- 
lation between gaugino masses. Mi : M2 : M3 = 
tti : a2 : as, is used as an additional assumption, 
where are the gauge coupling constants. All 
other SUSY masses which are needed for the cal- 
culations are treated as free parameters, limited 
from below by accelerator data. 

One can find the relevant calculations within 
this approach in refs. p3| , p^ and in the review p5[| 
(see also the references therein). There are large 
areas in neutralino parameter space where the 
neutralino relic density satisfies the relation (||). 
This is especially true for heavy neutralinos with 
> 100 - 1000 GeV, ref.[|6|. In these areas 
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there are good prospects for indirect detection of 
neutralinos, due to high energy neutrino radia- 
tion from Earth and Sun (see (2^:^^ and refer- 
ences therein) as well as due to production of an- 
tiprotons and positrons in our Galaxy. The direct 
detection of neutralinos is possible too, though 
in more restricted parameter space areas of light 
neutralinos (see review ]25|]). 

This model-independent approach is very inter- 
esting as an extreme case: in the absence of an ex- 
perimentally confirmed SUSY model it gives the 
results obtained within most general framework 
of supersymmetric theory. 

(ii) Strongly constrained models. This approach 
is based on the remarkable observation that in 
the minimal SUSY SU(5) model with fixed parti- 
cle content, the three running coupling constants 
meet at one point corresponding to the GUT mass 
Mgut- Because of the fixed particle content 
of the model, its predictions are rigid and they 
strongly restrict the neutralino parameter space. 
This is especially true for the limits due to pro- 
ton decay p K^v. As a result very little space 
is left for neutralino as DM particle. Normally 
neutralinos overdose the Universe (Jl^ > 1). The 
relic density decreases to the allowed values in 
very restricted areas where xX'^-nnihilation is ac- 
cidentally large (e.g. due to the exchange term 
- see ref. Thus, this approach looks rather 
pessimistic for neutralino as DM particle. 

In several recent works less restricted 

SUSY models were considered. In particular the 
limits due to proton decay were lifted. A GUT 
model was not specified or less restrictive S'O(IO) 
model was used Although, the neutralino 

can be heavy in these models, the prospects for 
indirect detection, including the detection of high 
energy neutrinos from the Sun and Earth, are 
rather pessimistic [Q. The direct detection is 
possible in many cases 34 , 3^] . 

(m) Relaxed restrictions. In section 5.3, follow- 
ing refs.|^^,|3^, we shall analyze the restrictions 
to neutralino as DM particle, imposed by basic 
properties of SUSY theory. As in many pre- 
vious works, a fundamental element of analysis 
is the radiatively induced EW symmetry break- 
ing (EWSB)lQ. However, some mass unifica- 
tion conditions at the GUT scale are relaxed (see 



ref.|36|). The powerful restriction from the no- 
finc-tuning condition is added. 

5.1. SUSY theoretical framework 

The basic element which should be used in 
the analysis is a supersymmetry breaking and in- 
duced by it (through radiative corrections) elec- 
troweak symmetry breaking |Q. We shall refer 
to this restriction as to the EWSB restriction. 

One starts with unbroken supersymmetric 
model described by some superpotential. It is 
assumed that local supersymmetry is broken by 
supergravity in the hidden sector, which commu- 
nicates with the visible sector only gravitation- 
ally. This symmetry breaking penetrates into the 
visible sector in the form of global supersymme- 
try breaking. More specifically it is assumed that 
the symmetry breaking terms in the visible sec- 
tor are the soft breaking terms given at the GUT 
scale ^ ^^gut by the following expression: 

Lsb = |0aP + mi/2^AQAa + 

a a 

+ Aniofy + BmofiHiH2 (11) 

where (j)a are scalar fields of the model (sfermions 
and two Higgses Hi and H2), Xa are gaug- 
ino fields, fy are trilinear Yukawa couplings of 
fermions and Higgses and the last term is an 
additional (relative to the superpotential term 
IJLH1H2) soft breaking mixing of two Higgses. 
Here and everywhere below we specify the scale at 
which an expression and parameters are defined. 

The soft breaking terms (|ll|) are described by 
5 free parameters : mQ,mi/2, B and /i. This 
implies the strong assumption that all scalars (pa 
and all gauginos Xa at the GUT scale have the 
common masses mg and ?tii/2, respectively. This 
assumption can be relaxed, as we shall discuss 
later. 

The soft breaking terms ( pT|) together with su- 
persymmetric mixing give the following potential 
defined at the EW scale at the tree level: 



V 



m2|i/i|2 + m2|i/i|2 



BfimQ{HiH2 + he) + 

d±A^\Hi\^^\H2\'r, 



(12) 
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The mass parameters mi and m2 at the GUT 
scale are equal to 



ml{GUT) = ml{GUT) = fi^ 



'0: 



(13) 



with fi defined at the GUT scale, too. The term fi^ 
in Eqs.(p^),(13) appears due to the mixing term 
I1H1H2 in the superpotential of unbroken SUSY. 

The radiative EWSB occurs due to evolution 
of TO_f/2j the mass of H2, which is connected with 
the upper components of the fermion doublets 
and in particular with t-quark. Because of the 
large mass of t-quark and consequently the large 



Yukawa coupling Yuh^ 



evolves from mi at 



the GUT scale to the negative value m'jj^ < at 
the EW scale. At this value the potential ( p^ ) 
acquires its minimum and the system undergoes 
the EW phase transition coming to the minimum 
of the potential. At EW scale in the tree approx- 
imation the conditions of the potential minimum 
(vanishing of the derivatives) give: 



tan^/? 



sin 2(3 ■ 



tan^ (3 - 1 



'Hi 



'H2 



2^2 



(14) 



With these equations we obtain one connec- 
tion between five free parameters describing the 
soft-breaking terms (pT|). Thus the number 
of independent parameters is reduced to four, 
e.g.mo,r7ii/2,v4, ^ (or tan/3). 

Using the renormalization group equations 
(RGE) one can follow the evolution of the scalar 
particles (Higgses and sfermions) and spin 1/2 
particles (gauginos) from the masses toq and 1111^2 
at the GUT scale to the masses at the EW 
scale. Analogously, the evolution of the coupling 
constants can be calculated. In particlular the 
masses of Higgses at EW scale are given by 



4. 



2 ' 6,™^ 



12^2 



c.A'm', + d,Amom,/2, 



(15) 



where ai,bi,Ci and di are numerical coefficients, 
which depend on tan (3. 

Equivalently, using Eqs.(O) one finds 



(16) 



where Ji are also numerical coefficients. 

Eq.(|l6|) allows to impose the no- fine-tuning 
condition in the neutralino parameter space. In- 
deed, one can keep large values of masses in the 
rhs of Eq([l^) only by the price of accidental com- 
pensation between the different terms. It is un- 
natural to expect an accidental compensation to 
a value less than 1% from the initial values. This 
is the no-fine-tuning condition. 
Naturally this condition is just the same as the 
one due to the radiative corrections to the Higgs 
mass. 

5.2. Restrictions: the price list 

Within the theoretical framework outlined 
above one can choose the restrictions from the 
following price list: 



• Soft breaking terms (11) and EWSB condi- 
tions (|lT 

• No-fine-tuning condition, 

• Particle phenomenology (constraints from 
accelerator experiments and the condition 
that the neutralino is the LSP), 

• Restrictions due to 6 ^ 57 decay, 

• Meeting of coupling constants at Mqut, 

• b — T and b — t — t unification, 

• Restrictions due to p — ^ Kv decay. 

Using some (or all) restrictions listed above one 
can start calculations for the neutralino as DM 
particle. The regions where the neutralinos are 
overproduced [fl^h^ > 1) must be excluded from 
consideration and the allowed region should be 
determined according to the chosen cosmological 
model (e.g. fl^h^ = 0.2 ± 0.1 for the CHDM 
model). For the allowed regions the signal for 
direct and indirect detection can be calculated. 

5.3. SUSY models with basic restrictions 

Accepting all restrictions listed above one ar- 
rives at a rigid SUSY model, with the neutralino 
parameter space being too strongly constrained. 
In ref.|^^ the SUSY models with basic restric- 
tions were considered. These restrictions are as 
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follows: 

(i)Radiative EWSB, (ii) No fine-tuning stronger 
than 1%, (iii) RGE and particle phenomenology 
(accelerator limits on the calculated masses and 
the condition that neutralino is LSP), (iv) Lim- 
its from b —>■ s"f decay taken with the uncertain- 
ties in the calculations of the decay rate and (v) 
0.01 < ^x^"^ < 1 as the allowed relic density for 
neutralinos. Rather strong restrictions are im- 
posed by the condition (ii); in particular it limits 
the mass of neutralino as < 200 GeV. 

At the same time some restrictions are lifted 
as being too model-dependent: (i) No restrictions 
are imposed due to p — s- Kv decay, (ii) Unifica- 
tion of coupling constants at the GUT point is 
allowed to be not exact (it is assumed that new 
very heavy particles can restore the unification), 



(iii) unification in the soft breaking terms (11) 
is relaxed. Following rcf.|3^ it is assumed that 
masses of Higgses at the GUT scale can deviate 
from the universal value mo as 



ml.(GUT)=ml{l + S,) (i = l,2). 



(17) 
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Figure 2. The neutralino parameter space for the 
mass-unification case di — 62 — and tan/3 = 8. 



This non-universality affects rather strongly 
the properties of neutralino as DM particle: 
the allowed parameter space regions become 
larger and neutralino is allowed to be Higgsino- 
dominated, which is favorable for detection. 

Some results obtained in ref. |Q are illustrated 
by Figs. 2-6. The regions allowed for the neu- 
tralino as CDM particle are shown everywhere by 
small boxes. 

In Fig. 2 the regions excluded by the LEP and 
CDF data are shown by dots and labelled as LEP. 
The regions labelled "fine tuning" have an acci- 
dental compensation stronger than 1% and thus 
are excluded. No-fine-tuning region inside the 
broken-line box corresponds to a neutralino mass 
< 200 GeV. The region "EWSB-fparticle 
phenom." is excluded by the EWSB condition 
combined with particle phenomenology (neu- 
tralino as LSP, limits on the masses of SUSY 
particles etc). In the region marked by rar- 
efied dotted lines neutralinos overdose the Uni- 
verse {il^h'^ > 1). The solid line corresponds to 
niQ = 0. The regions allowed for neutralino as 
CDM particle (0.01 < fl^h'^ < 1) arc shown by 



small boxes. As one can see in most regions the 
neutralinos are overproduced. The allowed re- 
gions correspond to large XX annihilation cross- 
section (e.g. due to Z°-pole). 
Fig. 3 and Fig. 2 differ only by universality: in 
Fig. 3 (5i = (52 = (mass-unification), while in 
Fig. A61 = —0.2 and 62 = 0.4. The allowed region 
in Fig. 3 becomes much larger and is shifted into 
the Higgsino dominated region. Figs. 4 and 5 are 
given for tanf3 = 53. This large value of tan/3 
correspond to b — t — t unification of the Yukawa 
coupling constants. Again one can notice that 
in the mass-unification case (^i = = 0) only 
a small area is allowed for neutralino as CDM 
particle, while in non-universal case (5i — 0,62 = 
—0.2) the allowed area becomes larger and shifts 
into the gaugino dominated region. 

In Fig. 6 the scatter plot for the rate of direct 
detection with the Ge detector |^ is given for 
the non-universal case {61 — 0,62 ^ —0.2) and 
tan/3 = 53. We notice that, for some configura- 
tions, the experimental sensitivity is already at 
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Figure 3. Case 5i = -0.2, 82 = 0.4 and tan /? = 8. 



Figure 4. Mass-unification case and tan /3 = 53 



the level of the predicted rate. 

6. Conclusions 

1. The density of CDM needed for most cosmo- 
logical models is given by ilcDAdh^ = 0.2 ± 0.1. 
There are four candidates for CDM, best moti- 
vated from point of view of elementary particle 
physics: neutralino, axion, axino and majoron. 

2. There are different approaches to study the 
neutralino as DM particle in SUSY models with 
R-parity conservation. 

In the phenomenological approach, apart from 
the LEP-CDF limits, very few other constraints 
are imposed. In the Minimal Supersymmetric 
Model only one GUT relation between gauginos 
masses, Mi : M2 : M3 = ai : : as, is used. 
While the coupling constants are known, the mass 
parameters needed for calculations are taken as 
the free parameters. Many allowed configurations 
in the parameter space give the neutralino as DM 
particle with observable signals for direct and in- 
direct detection (antiprotons and positrons in our 



Galaxy and high energy neutrinos from the Sun 
and Earth). 

The other extreme case, the complete SUSY 
SU(5) model with fixed particle content, with 
meeting of coupling constants at the GUT point 
and with the constraints due to proton decay, 
leaves very little space for the neutralino as DM 
particle. 

In the third option the SUSY soft breaking 
terms (|l^) and induced EW symmetry breaking 
are used as the general theoretical framework. 
Combined with a no-fine-tuning condition this 
framework already gives essential restrictions. In 
particular, fine tuning allowed at the level larger 
than 1% results in the neutralino being lighter 
than 200 GeV. Here the neutralino is gaugino 
dominated, which is unfavorable for direct detec- 
tion. If we employ further other restrictions, such 
as exact unification of gauge coupling constants 
and soft-breaking parameters at the GUT scale, 
b —f limit and b — t unification, the model 
becomes as rigid as the one considered above. 

If, on the other hand, we relax some con- 
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Figure 5. Case Si — 0, S2 — —0.2 and tan/3 = 53. 



Figure 6. Rate for direct detection {Si — 0, S2 
-0.2 and tan/3 = 53). 



straints, 6.17. by assuming non- universality of the 
scalar mass term in Eg. (pi]), then even in the case 
of 1% no-fine-tuning condition the neutralino can 
be the DM particle in a large area of the param- 
eter space and can be detected in some parts of 
this area in direct and indirect experiments. 

The neutralino is an observable particle. It can 
be observed directly in the underground experi- 
ments, or indirectly, mostly due to the products 
of neutralino-neutralino annihilation. 

3. In the framework of supersymmetric the- 
ory, the PQ-mechanism for solving the problem of 
strong CP violation, results in a Goldstone super- 
multiplet which contains axion and its fermionic 
superpartner, axino. Both of them can be CDM 
particles. The most important parameter here is 
the scale of PQ symmetry breaking fpq which 
is observationally constrained as 2 • 10^ GeV < 
fpQ < 8 • 10" GeV. The axion can be the CDM 
particle, if its mass is 10~^ — 10~^ eV ( the corre- 
sponding values of fpg are between 2 • 10^ GeV 
and 6 • 10^" GeV). For larger values of fpq the 
axino can be CDM particle if a mechanism of 



SUSY breaking provides its mass within the in- 
terval 10 — 100 keV. The axino can provide both 
CDM and HDM components needed to fit the cos- 
mological observations. The axion can be directly 
observed (e.g. in microwave cavity experiments) 
while the axino dark matter is practically unob- 
servable. 

4. The Majoron with keV mass can be the 
warm DM particle which explains the galactic 
scale (M ~ in the structure formation 

problem. The decay of the Majoron to two pho- 
tons can produce an observable X-ray line in the 
cosmic background radiation. 
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